#install.packages("ggseqlogo")
#install.packages("cowplot")
library(ggseqlogo)
library(stringr)
library(ggsci)
library(tidyverse)
library(ggprism)
library(cowplot)

setwd("C:\\Users\\dell\\Desktop\\")

seq_dna <- read.table('seqlogo.txt', header = T)

# method = 'bits'
p1 <- ggseqlogo(as.character(seq_dna$uORF.T2), method = 'prob',
                font="roboto_slab_light", col_scheme=mycol1) +
  theme_bw() +
  scale_y_continuous(labels = scales::percent)
p1

p2 <- ggseqlogo(as.character(seq_dna$uORF.nT2), method = 'prob',
                font="roboto_slab_light", col_scheme=mycol1) +
  theme_bw() +
  scale_y_continuous(labels = scales::percent)
p2

plot_grid(p1, p2, ncol = 1, align = 'v')
# Create custom colour scheme
cs2 = make_col_scheme(chars=c('A', 'T', 'C', 'G'), values=1:4)

# Generate sequence logo
ggseqlogo(seq_dna$uORF.T, col_scheme=cs2, method = 'prob')


### 

mycol1 <- make_col_scheme(chars=c('A', 'T', 'C', 'G'),
                          cols=c('#5c8d89', '#74b49b', '#ffbcbc', '#a7d7c5'))

ggseqlogo(as.character(seq_dna$uORF.T), method='prob', col_scheme=mycol1)







